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Abstract 

We discuss two polynomial bi-Hamiltonian structures for the generalized integrable 
Chaplygin system on the sphere S 2 with an additional integral of fourth order in momenta. 
An explicit procedure to find the variables of separation, the separation relations and the 
transformation of the corresponding algebraic curves of genus two is considered in detail. 
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1 Introduction 

We address the problem of the separation of variables for the Hamilton- Jacobi equation within 
the theoretical scheme of bi-hamiltonian geometry. The main aim is the construction of different 
variables of separation for the given integrable system without any additional information (Lax 
matrices, r-matrices, links with soliton equations etc.) 

Different variables of separation may be useful in different perturbation theories, as well as 
distinct procedures of quantization and various methods of qualitative analysis etc. From math- 
ematical point of view the notion of different variables of separation allows us to study relations 
between distinct algebraic curves associated with the corresponding separated equations. Such 
relations give us a lot of new examples of reductions of Abelian integrals and, therefore, they 
may be a source of new ideas in number theory, algebraic geometry and modern cryptography. 

The paper is organized as follows. In Section II, the necessary aspects of bi-hamiltonian 
geometry are briefly reviewed. Then, we discuss a possible application of these methods to 
calculation of the polynomial bi-hamiltonian structures for the given Chaplygin system. In 
Section III, the problem of finding variables of separation and corresponding separation rela- 
tions is treated and solved. Some applications of real and complex variables of separation are 
discussed in the final Section. 

2 The bi-hamiltonian structure of the Chaplygin system 

In order to get variables of separation according to general usage of bi-hamiltonian geometry 
firstly we have to calculate the bi-hamiltonian structure for the given integrable system with 
integrals of motion Hi, ... , H n on the Poisson manifold M with the kinematic Poisson bivector 
P and the Casimir functions C x , . . . , C k [S E QJ OS E] ■ 

Let us consider the generalized Chaplygin system jH 03 ED with the following integrals of 
motion 

H x = Jl + Jl + 2Jl + c 2 {x\-xl) + ^, c 2 ,c 4 eK, (2.1) 

x 3 

H 2 = (ji+Ji + %) +2c 2 xl(Jl -J%) + <%b\. 
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In this standard coordinates J = (Ji, J 2 , J3) and x = {x\, x 2 , X3) on the Euclidean algebra e*(3) 
the Poisson bivector looks like the following antisymmetric matrix 
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It has two Casimir elements 



PdC h2 = 0, 



Ci 



k=l 



Co 



(x, J) = ^X k J k 



(2.2) 



fc=i 



At C 2 = {x, J) = this integrals of motion are in involution 

{H U H 2 } = {PdH u dH 2 )=0. 

and the corresponding symplectic leaves are equivalent to cotangent bundle T*S 2 of the sphere 
with the radius a = \x\. We will consider such symplectic leaves only, and, therefore, all the 
formulae below hold true up to C 2 = 0. 

Following to [TH [T31 HH E] we will suppose that the desired second Poisson bivector is the 
Lie derivative of P along some unknown Liouville vector field X 



which has to satisfy to the equations 

{H u H 2 y ={P'dH!,dH 2 ) 



P' = C X {P) (2.3) 



0, [P , ,P'] = {C x (P),Cx(P)]=0, (2.4) 



where [., .] means the Schouten bracket. 

Obviously enough, in their full generality equations (I2.4[) are too difficult to be solved 
because it has infinitely many solutions 114) . In order to get some particular solutions we 
will use the additional assumption 

P'dCx.2 = 0, (2.5) 

and polynomial in momenta J ansdtze for the components X J , j = 1, . . . ,6, of the Liouville 
vector field X = £ X* d 3 ; [HI H2 [13] . 

Substituting this ansdtze into the equations fj2.4H2.5f) and demanding that all the coefficients 
at powers of J vanish one gets the over determined system of algebro-differential equations on 
functions of x which can be easily solved in the modern computer algebra systems. 

For the linear ansatze one gets the trivial solutions P 1 = P only. The quadratic ansatze 
yields first non-trivial solution 



X = 



w 3 

a 2 x 3 



(wi,W2,w 3 , 0,0,0) - ^(0,0,0,^3,0,-^1), 



W = J x x, 



where W is a cross product of J and x. This Liuville vector field gives rise to the following real 
Poisson bivector 



P' 
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X2J3 

xi. J '3 
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(2.6) 



2 



For the cubic ansatze we get one real Poisson bivector compatible with (|2.6|) and two new 
hermitian conjugated bivectors P' and P'*. This new and more cumbersome solution reads as 

P[ 2 = x 3 {W 3 + ix 3 J 3 ), P[ 3 = ~x 2 (W 3 + ix 3 J 3 ), J* 3 = xi(W 3 + ix 3 J 3 ) 

PL = J 2 w, - ix 3 j x j 3 - ^t 2 x 3 {x 2 j 3 - wo + xl{c2X \ - ca) 

x 3 

p{ 5 = -(* + i W - iftf - 4) + W»M ~ Wi) X3 - " 2x ^ x * + c *> 

P{ e = -iJsWr + " 2(r - iC2( f - 4X1X2)) - - «» ~ h») + ^ 

P^ 4 = — i( Ji + iJ 2 )W 2 - f ( J? - Jf ) + ^(xi J, + W 2 ) - ^ + 2 ^1^4 ~ 



L 3 

„2l„ „4 



^25 = -JiWa - i Xl J 2 J 3 - iV&x 3 ( Xl J 3 + W 2 ) - Xx{ - C2X \ + C4) (2.7) 

x 3 

v> ■7TI/ ixi(T-ic 2 (ia;| +4xix 2 ) v ic 4 3Ji 
P 26 = -1J3W2 jr- 2 V^(a;iJ3 + W2){xi - 1X2) — 

Z x 3 

PL - J*w 3 + iX2(J [ + J|) + ^ 3 ws + iX2{C2 **f C4) 

Pgg = -U3W3 - y/c2~W 3 (xi - IX 2 ) - 1C 2 XiX 2 X 3 



. i(J? + Jf) J 3 , — / 2 , . . (xi — ix 2 )c 4 T/T . T , 
^45 = caa§(gi+ia;2) + - ^ x 3 J 3 {J x -U 2 ) 



\xi{Ji + 2iJ 2 )(c 2 a;3 + c 4 ) ix 2 (J 2 + 2iJ0( c 2^3 - c 4 ) 

^46 = - ^^^(^i + i^ 2 ) - (si - '^ 2 )(^J 2 J 3 - ~~^r) 
a;ix 2 ( 2c 2 J 2 + —4— ) — (a; 3 + zx 2 ) I c 2 - — r 



n/ iJiT — / . . . . / C 4 Xl 

^56 = q h VC2 |cixia;3(a;i + ix 2 ) - {xi ~ ix 2 )[ J\ J 3 3- 

z \ \ x 3 

+ x lX2 [2c 2 J 1 - 1 -^-) - 1 A(xl + 2x 2 1 ) (c 2 + ^ 



.. x 3 

Here T = Jj 2 + Jf + 2J| and i = \/— I. This cubic Poisson structure may be rewritten in 
lucid form by using 2x2 Lax matrices for the Chaplygin system [5j El] and the bi-hamiltonian 
structure associated with the reflection equation algebra [15] . 

The quartic ansatze yields a lot of solutions, which will be classified and studied at a future 
date. 

To sum up, using applicable ansatze for the Liouville vector field X we get a real relatively 
simple quadratic bivector (|2.6[) and a more complicated complex cubic bivector (|2.7p . Modern 
computer software allows to do it on a personal computer wasting only few seconds. The 
application of this Poisson bivectors will be given in the next section. 
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3 Variables of separation and separation relations 



The second step in the bi-hamiltonian method of separation of variables is calculation of canon- 
ical variables of separation (qi, . . . , q n ,pi, . . . ,p n ) and separation relations of the form 



4>i(qi,Pi,Hi,...,H n ) = , i= l,...,n, with det 



^0. (3.1) 



The reason for this definition is that the stationary Hamilton- Jacobi equations for the Hamil- 
tonians H t can be collectively solved by the additively separated complete integral 

n 

W(q x ,..., )=^W i {q i ;a l ,...,a n ), (3.2) 

i=l 

where Wi are found by quadratures as solutions of ordinary differential equations. 

According to [31 [T21 [T7] , separated coordinates qj are the eigenvalues of the control matrix 
F defined by 

F'dH = P(FdH). 

Its eigenvalues coincide with the Darboux-Nijenhuis coordinates (eigenvalues of the recursion 
operator) on the corresponding symplectic leaves. Using control matrix F we can avoid the 
procedure of restriction of the bivectors P and P' on symplectic leaves, that is a necessary 
intermediate calculation for the construction of the recursion operator [3J . 

Recall that independent integrals of motion (Hi, . . . , H n ) are Stackel separable if the cor- 
responding separation relations are given by the affine equations in Hj, that is, 

n 

}^ Sjj(pi, qi)Hj - Uj(pi, qi ) = , i=l,...,n, (3.3) 

3=1 

with S being an invertible matrix. Functions Sij and Ui depend only on one pair (pi,qi) of 
canonical variables of separation, thus it means that 

{Sik,qj} = {Sik,Pj} = {Sik,Sj m } = 0, i^j, (3.4) 

and 

{Ui, qj } - {Ui, Pj } = {Ui, Uj} = 0, i ? j. (3.5) 

In this case S is called a Stackel matrix, and U - a Stackel potential. 

For Stackel separable systems the suitable normalized left eigenvectors of control matrix 
F form the Stackel matrix S [3] 

F = S- 1 diag (qi,...,q n )S. 

However, we have to notice that the definition of Stackel separability depends on the choice 
of H. Indeed, if (Hi, . . . ,H n ) are Stackel-separable, then Hi = H(Hi, . . . , H n ) will not, in 
general, fulfill relations of the form (|3.3p . 

So, the main problems are the finding of the conjugated momenta pi, such that {pi, qj} = 5ij 
and the construction of separation relations (f>j p.l[) for generic non-Stackel separable systems. 
We show below how we can solve these problems using the same control matrix F with the 
addition of some useful observations. 

3.1 The quadratic Poisson bivector 

According to [3] , the bi-involutivity of integrals of motion 

{Hi,H 2 } = {H 1 ,H 2 }' = 
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is equivalent to the existence of the non-degenerate control matrix F such that is 

2 

P'dH = P(PdH) , or P'dH, = P ^ dHj, i=l, 2. 
For the quadratic in momenta Poisson bivector P' (I2.6P control matrix F reads as 



(3.6) 



F = 



( 



\ 



1 f J 2 + Jf 

2 



c 4 

C2 - — 

2=3 



1 



2^-71) + ^!-^ ~ 



1 /J 2 + J; 



C4 

C 2 + — 

X t 



(3.7) 



The eigenvalues of this matrix are the required variables of separation gi j2 



det(P - XI) = (A — <?i)(A — q 2 ) = A 2 



c 2 A + 



c 2 jf 



(3.8) 



If the corresponding separated relations are affine equations in H\2 then the suitable normalized 
left eigenvectors of F form the Stackel matrix S 

F = S~ 1 diag(q 1 ,q 2 )S. 

In our case the matrix of normalized eigenvectors 

2c 4 



S 



si s 2 

1 1 



Sl,2 



2 ± 2 J (J* + 4f + 2c 2 xl(Jl - Jf) + c\x\ 



(3.9) 



does not form the standard Stackel matrix because 

{QhSj} 7^0, and {si,s 2 }^0. 

It means that the underlying separation relations do not form the Stackel affine equations (|3 . 3[) 
in Pi j2 . Substituting functions on the integrals of motion Pi j2 = fi,i{Hi, H2) into the equation 
(EH) ' 

P'dH = P(PdH) , 

we can try to get a new control matrix F which satisfies the Stackel properties (13.4[) . In the 
Chaplygin case it is a very simple calculation which yields to the following results 



H 1 =H U 



Ho — Ho 



and 



F21 



F\2 — 



2C 2 J( - J; 



2 (^l ^X*^! + ^2 + ^ Jf 



c 2 ( x l x 2 + x 3)( x l x 2 X V) 



F22 



Jf c 2 (x 2 -x 2 



so that 



5 = 



Sl s 2 

1 1 



Si, 2 = 2[H 1 - 



C4 



± 2^ + Jf)3 + 2c 2 ^(J 2 - J 2 ) + c\x% . 



5 



The Stackel conditions (|3.4|) are fulfilled and, therefore, Si is a function on the separated 
coordinate q\ defined by (|3.8p and the unknown conjugated momenta p\. It is easy to see that 
the recurrence chain 

4>i={si,qi}, 02 = {<pi,qi }, 4>i = {<j>i-i,qi} (3.10) 

breaks down on the third step 03 = 0. It means that Si is the second order polynomial in the 
momenta p\ and, therefore, we can define this unknown momenta in the following way 



k _ _X3_ (x 2 J 1 - x 1 J 2 )qi + c 2 x x J 2 
b 2 ~ 2 (J 2 + Jf )<?! - c 2 J 2 



YL ^3 \^ 2 ^i 11 T /oil's 

= — = ~ — 7 72 T m7 7^2 ( 3 - n ) 



up to the canonical transformations pi — > pi + g(qi)- 

The similar calculation for S2 yields definition of the second momenta 



£3 (x 2 Ji - xiJ 2 )q 2 + c 2 xiJ 2 

P2 ~~Y (J* + Jl) q2 -c 2 Jl ■ (3 - 12) 



So, one gets canonical transformation from initial physical variables x, J to the variables of 
separation p, q (|3.8I3.11I3.18[) . 

The inverse mapping looks like 



<7i<72 T _ jqm 7 _ 2^Jq 1 q 2 {qi - c 2 )(c 2 - g^jjpi - p 2 ) 



J\ = \ qi + qi - c 2 x 3 , J 2 = J — x 3 , Jz = 



where x% = ^/ a 2 — x\ — x\ and 



C2 V c 2 91 - 92 

(3.13) 



., /<7i<72 Pi9i - P292 - c 2 (pi - p 2 ) (qi ~ c 2 ){c 2 - q 2 ) piqi - p 2 q 2 

xi=2* , x 2 — -2\ 



c 2 91-92 V ° 2 qi ~ 92 

Using the elements of the Stackel matrix 

Si,2 = 2(8(c 2 - qx <2 ) gi, 2 Pi, 2 ~ a 2 ( c 2 - 2q 1 ^ 2 )j 
we can easily derive the required affine Stackel separated relations 

S 2 

s k H x + H 2 = ^-- 2c 4 (c 2 - 2q k ) , k = 1, 2. 

If we come back to initial integrals of motion Hi 2 then these separation relations go over 
to the equation 

C: $(q,p) = (8q(c 2 -q)p 2 -a 2 (c 2 -2q)-H 1 + ^H 2 ~^ (3.14) 
x (8q(c 2 - q)p 2 - a 2 (c 2 - 2q) -E x - JW 2 \ - 2c 4 (c 2 - 2q) = . 



Proposition 1 XTie variables of separation {qi,Pi) He on the hyperelliptic algebraic curve C of 
the genus two g = 2 defined by \3. 11$ and the equations of motion are linearized on its Jacobian. 

Remark 1 The separation relations (|3.14p may be obtained in a framework of the Stackel 
formalizm (|3.3|) by using generalized Stackel matrix S (|3.9|) . whose entries Sij depend on {qi,Pi) 
and on the integrals of motion 

Si = Si(pi,qi) - 2Hi 

Remark 2 At C4 = we reproduced the Chaplygin result so that (qi,Pi) lie on a pair of the 
elliptic curves of genus one 

Ci, 2 : 8q i c 2 -q)p 2 -a 2 (c 2 -2q)-H 1 ±^H 2 ~ = 0. (3.15) 

Remind that every algebraic curve of genus one is isomorphic to a real torus. 
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3.2 The cubic Poisson bivector 

For the cubic in momenta Poisson bivector P' (|2.7[) the entries of control matrix F look like 

Fn = ~iT + ^/c 2 '(x 3 (J 1 ~U 2 )~2(x 1 -ix 2 )J 3 ) +2c 2 x 1 x 2 - 1 -^-, 

x 3 

F 21 = -i (.Jf + J 2 + -^0 (jl + jf + f| + 2^x 3 i( Ji - iJ 2 )^ + ic 2 ^ 
+ 2c 2 ^(^.T3(Ji + iJ 2 ) + 2Ji J 2 ) , 

F\ 2 = — , i 7 22 = . 

The separated variables Ai j2 are the roots of the characteristic polynomial 

det(F - XI) = (A - Ai)(A - A 2 ) = A 2 - F M A - ^F 21 . (3.16) 

In the next step we have to find the conjugated momenta to these Darboux-Nijenhuis coordi- 
nates. 

In [5|ll0j we've found the separated coordinates Qi, 2 and the corresponding momenta Pi i2 
for the Kowalevski-Goryachev-Chaplygin gyrostat using 2x2 Lax matrix, its Baker- Akhiezer 
vector-function and the reflection equation algebra. In the framework of Sklyanin formalism 
the separated coordinates are the poles of Baker- Akhiezer function with standard simplest 
normalization, whereas the separated momenta are expressed through the values of the elements 
of the Lax matrix in the poles. 

It is easy to see that the Darboux-Nijenhuis variables Ai. 2 (|3.16p are related with the poles 
<3i, 2 of the Baker- Akhiezer function by the following point transformation 

Ai,2 = ~Q 2 , 2 . (3.17) 

In this case the additional knowledge of the Lax matrix allows us to introduce conjugated to 
Afe momenta 

^ k 21 u=Q k 

where B(u) is the diagonal element of the corresponding Lax matrix |10| 



(3.18) 



B(u) = ~(xi - ix 2 )u 3 + (j2J 3 - iy/c2~(xi - ix 2 )){xi - ix 2 ) - (Ji - iJ 2 )x 3 ) 
2x 3 J 3 (J\ — iJ 2 ) + ( J\ + Jf + -| ) (xi — ix 2 ) — c 2 (xi + ix 2 )xl ) u 



v 2 



ic 3 2 /2 xl + c 2 xl( Ji + iJ 2 ) + iyfcxl{Jx - iJ 2 ) 2 + x 3 ( Jl + J| + % ) (Ji - iJ 2 ) 
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Moreover, it allows us to prove that variables of separation (Xi,[Ai) lie on the Jacobian of the 
other hyperelliptic curve of genus two 

C: <f>(z,X) = z 2 + (4X 2 +4LffiA- H 2 )z + (2a 2 A + i c 4 ) 2 c 2 = 0, (3.19) 



where we put z = — i v / c 2 e 2/iV ^\ 

So, we know the answer and, therefore, we could try to guess how to get this information 
without using our knowledge of the Lax matrix. Namely, the normalized left eigenvectors of F 
form the Stackel matrix 



P = ^ 1 diag(A 1 ,A 2 )^, S=(-f' -f 2 
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In order to get conjugated to Ai, 2 momenta mi,2 we could use the Stackel potentials (|3.5p 

U h2 = -4iAi, 2 ifi + K 2 , such that {{/i, J7 2 } = 0, {U h Xj} = 0, i^j. 

As it has been stated above we could study the following recurrence chain of the Poisson 
brackets 

<h. = {XiM, 2 = {Ai,0i} > ..., (f>i = {Xi,<f>i-\}. (3.20) 
This chain is a quasi-periodic chain 

03 = 8iAi0i. 

It means that the Stackel potential U\ is some trigonometric function on momenta Mi and, 
therefore, we could define this desired momenta in the following way 

Mi = /(Ai)ln(V8iAi>i + <£ 2 ) 

up to canonical transformations /ii — > Ml + g(Xi). Here function f{X\) is determined from the 
canonicity of the bracket {Ai, Mi} = 1. 

Summing up, the central idea of the proposed construction is the observation of quasi- 
periodicity or of the break of the recurrence chain of the Poisson brackets (|3.10l) or (I3.20j) . 
We have to point out again that all the appropriate tedious calculations may be done on the 
personal computer in a few seconds. So, it is a real way to get the variables of separation and 
the separated relations for the given integrable system. 

4 The quadratures 

According to the Liouville theorem the existence of n independent integrals of motion Hi in 
the involution {Hi, Hj} — guarantees that equations of motion may be solved in quadratures. 
Usually we fit some additional requirements on these quadratures, as for example, the solutions 
of equations of motion have to be single valued real functions on the real time variable t. It is 
necessary for the qualitative analysis of motion, for topological analysis, for perturbation theory 
etc. 

Of course, we prefer to have real solutions represented by an explicit closed-form expression 
right away. However, in nine times out of ten in order to get such single- valued solutions we 
have to start with the solution of the Jacobi inversion problem. For example, for the geodesic 
motion on an ellipsoid and for the Kowalevski top we have to find variables of separation si, 2 (i) 
from the equations 

ds r S2 ds 

sds f S2 sds 

where P is a polynomial of degree 5 or 6, and /3i,2 are two constants of integration. Secondly, 
we have to express the initial real variables via variables of separation si j2 (t). 

In the previous section we get the real variables of separation (q,p) and the complex vari- 
ables of separation (A, /i) or (Q, P). These variables are related by the canonical transformation 

Ai,2 = Xi :2 (qi,q2,Pi,P2), Hi,2 = Mi, 2(qi, q2, Pi, P2) 

which may be obtained directly from the definitions (|3.16I3.18|) and the mapping (I3.13[) . Using 
the separation relations p,14[) we can rewrite this transformation as a quasi-point canonical 
transformation [B] 

Ai,2 = Ai,2(9i) Q.2) Hi,H 2 ) 
which relates the Jacobian of the hyperelliptic curve C f|3 . 14[) with the Jacobian of the hyperel- 
liptic curve C p,19[) . However, we can not rewrite it as a rational mapping ip : (p, q) — > (A, p,), 
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i.e. it does not be well-studied cover of the hyperelliptic curve of genus two even at C4 = 0, 
see [1] and references within. We suppose that these Jacobians are non-isogeneous in Richelot 
sense [5] as well. Any further inquiry of this relation goes beyond the scope of this paper. 

We have no right to escape complex variables and prefer the real variables of separation 
only. It is how these variables are employed that determines the good or evil. As often as 
not the equations of motion are linearized on the abelian variety, which is roughly spiking 
the complexified of the corresponding Liouville real torus, see more detailed discussion in [2]. 
Moreover, we remind that Lyapunov could improve the Kowalcvski result and solve the problem 
pointed out by Painleve by using the complex time variable t [7]. The Ziglin method [19], 
the theory of algebraically integrable systems [16] and some other modern theories deal with 
the complex analytical Hamiltonian systems only. Quantum mechanics is formulated over the 
complex field as well. 

4.1 The real variables of separation 

In this case equations of motion reads as 

dqi „ , . (. c 4 (gi - gg) 



dt ' \ (4gi(c 2 - qi)p\ - 4q 2 (c 2 - q 2 )p\ - (q% - <Z2)a 2 )/ ' 

(4.21) 

d q2 i s , c 4 {qi - q 2 ) 

- 12{C 2 - q 2 )p 2 1 + 



dt \ (4qi(c 2 - qi)pj - Aq 2 (c 2 - q 2 )p 2 - (qi - q 2 )a 2 ) 

At C4 = we have an accidental degeneracy of the genus two algebraic curve (|3.14j) to a product 
of two elliptic curves (I3.15[) . Namely, if C4 = then 

dqi *» 2dt 



8<7i(c 2 - q\)pi 8q 2 (c 2 -q 2 )p 2 

and we have to find functions qi j2 (t, ct\ t2 , fi\^ 2 ) from the two independent equations 

91 dq , f q2 dq 

, =Px+t, and / - =P 2 -t. 

Here we fix the values of the integrals of motion H\ 2 — a\^ 2 , and 

Pfc(g) = 8gi >2 (c 2 - qi,2)(a 2 (c 2 - 2q ia ) + a x =p 

are the cubic polynomials defining two elliptic curves (|3.15p . According to pQ in this case we can 
get solutions of Xk (t) and Jk (t) in an explicit and closed form using the Weierstrass p function 
and its derivative. 

At C4 ^ we have to solve the standard system of the Abel-Jacobi equations 

d<1 ■ ^ ^ fc + t, 



and 



91 



16g(c2-g)P(«) l6q(c 2 -q)P{q) 
dq 

00 32q{c 2 ~ q) (a 2 (c 2 - 2q) + ai - 8<?(c 2 - q)P 2 (q)) P(g) 

d A = p 2 

00 ?,2q{c 2 ~q)(a 2 {c 2 -2q)+a l -8q(c 2 -q)P 2 {q)Y{q) 

where P(q) means the solution of the equation p) = (|3.14[) with respect to p. For example, 
we can solve this equations numerically by using the Richelot approach [S] . 
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The second part of the Jacobi separation of variables method consists of the construction 
of the integrable system starting with some known separated variables and some arbitrary 
separated relations [T7]- As an example, if we consider new separation relations 

$(p,g) = <S>{p,q)~c 5 q 2 =0, 

where <J>(p, q) is given by (|3.14[) , then we get the following generalization of the Chaplygin 
system 

Hi = H x - + J 2 + c 2 x\) . 

We could not get any physically interesting systems by using these real variables of separation. 

4.2 The complex variables of separation 

In this section we will work with variables {P,Q) instead of (ji,X), i.e. we will consider a 
ramified two-sheeted covering given by point transformation p. 171) . 

The poles of the Baker- Akhiezer function lie on the Jacobian of the 2x2 Lax matrix 
spectral curve p.0.. This curve is defined by the equation with the real coefficients 

u) = y 2 ~ (u 4 ~ 2H lU 2 + H 2 )y + c 2 (a 2 u 2 - c 4 ) 2 - 0. (4.22) 

However, in order to get the separation relations we have to substitute the complex functions 
Qi,2 and Pi,2 on the initial real variables into this equation 



w = Qi,2, y = c 

The Abel- Jacobi equations have the standard form 



- 2iPl,2 



fQi rQ? rQi 

t + fa.= I 0,1+ I fil, 02= / ft 2 +/ ^2 



OC 



where 



d<f>(y,u)/dy d<P(y,u)/dy 

After the solution of the Abel- Jacobi equations we have to express the real initial variables x, J 
in terms of this complex solutions Qi, 2 (t) and Pi, 2 (t). 

On the other hand, these complex variables of separation are very useful for the construction 
of other integrable systems. Namely, substituting Qi t2 and Pi y2 into other separation relations 
we can get the Hamilton function for the Kowalevski-Goryachev- Chaplygin gyrostat [SJ [TU] 

Hi = Jl + Jf + 2J| + pJ 3 + cixi + c 2 {x\ - x\) + c 3 xia;2 + % (4.23) 



x 



3 



after some additional canonical transformations. Moreover, using these complex variables we 
can easily get quantum counterpart of the Chaplygin system and its generalizations [5]. 



5 Conclusion 

Starting with integrals of motion for the generalized Chaplygin system we have found the 
two polynomial in momenta Poisson bivectors (|2.6j) and (|2.7p . which are compatible with the 
canonical Poisson bivector on cotangent bundle T*S 2 of the two-dimensional sphere. 

An application of the corresponding control matrices allows us to get two families variables 
of separation and of separated relations using methods of bi-hamiltonian geometry only. The 
solutions of equations of motion in these variables are briefly discussed. 

The proposed approach may be useful for investigations of other integrable systems on the 
sphere with integrals of motion higher order in momenta |18j . for instance, the search another 
real variables for the Kowalevski-Goryachev- Chaplygin gyrostat (|4.23[) . 
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